Thermodynamics of a mixed quantum-classical Heisenberg model in two dimensions 
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, We study the planar antiferromagnetic Heisenberg model on a decorated hexagonal lattice, 

involving both classical spins (occupying the vertices) and quantum spins (occupying the middle of 
the links) . This study is motivated by the description of a recently synthesized molecular magnetic 
, compound. First, we trace out the spin | degrees of freedom to obtain a fully classical model with 

an effective ferromagnetic interaction. Then, using high temperature expansions and Monte Carlo 
simulations, we analyse its thermal and magnetic properties. We show that it provides a good 
quantitative description of the magnetic susceptibility of the molecular magnet in its paramagnetic 
phase. 

The Heisenberg model jlj has a long history and has been extensively studied throughout these last thirty years. 
While it is exactly solvable in one dimension B in some of its versions, it is only through approximate methods that 
quantitative information can be obtained in higher dimensions. High and low temperature expansions I^Q, Monte 
Carlo simulations and renormalisation group calculations have been widely developed and give now a precise 
account of the critical regime of the model. However, little has been done in the various specific contexts which are 
£f) , now realized in the magnetic molecular materials. 

For instance, the compound (NBu4)2Mn 2 [Cu(opba)]3-6DMSO-H 2 0, recently synthesized by H.O. Stumpf et al [||, 
exhibits a transition at T c — 15K towards an ordered state. The structure of this material can be schematically 
described by a superposition of layers of hexagonal lattices with the Mn 77 ions occupying the vertices and the Cu 77 
ions occupying the middle of the links, as shown in Fig. 1. The interplane-coupling is small, so that the spin system can 
be considered two-dimensional. In the plane, the nearest neighbour Mn-Cu ions interact through an antiferromagnetic 
coupling. It is interesting to determine the extent to which such a simple microscopic model with no other interaction 
. included, can quantitatively describe the magnetic and thermal properties of such a complex molecular architecture. 
Of course, the isotropic 0(3) model is critical only at zero temperature [jlO| and the symmetry breaking at T c = 15K 
has presumably its origin in a slight spin anisotropy and/or a small interplane coupling. However, one expects for 
O ■ T>T C that the properties of the material are well described by the two-dimensional isotropic antiferromagnetic spin 
| ^ - spin | interaction. This is the problem we investigate in this paper. 

We denote by S^ Mn ^ the spin | operator associated with the Mn ion at site j, and by s| Cu ^ the spin | operator 
corresponding to the Cu ion at site i in the middle of a link of the honeycomb lattice. The antiferromagnetic interaction 



' is represented by the Heisenberg hamiltonian 
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n = jJ2 S * (CU) ' S f n) - H £ S; (Mn) - g^ B H £ St (Cu) (1) 

<i,j> j=l i=l 

where J is positive, H is the external magnetic field, < i,j > stands for a pair of nearest neighbour spins, Ns is the 
number of sites and Nl is the number of links on the honeycomb lattice (Nl = 3/2Ns). The spin | operator can 

be approximated by a classical spin S.s where s is a unit classical vector and S = y§(§ + 1), whereas the spin | 

operators are expressed in terms of the Pauli matrices, S^ Cu ^ = | tr. Since the quantum spin sites are not directly 
coupled to each other, one can trace out the quantum spin dependence to get a completely classical partition function 

Z(T,H) = Jvs | Y[ 2cosh(||-i/3J5(si + s/) +a 2 He z \\)\ exp ^ai-ff ^ sfj (2) 

where we have defined ai = S f3gxfj,B, ct2 = \^Qi^b^ *Ds = rij=i sm ^j d6j d<pj and ||X|| stands for the length of 
vector X. The indices i and j now label the classical spins located at the vertices of the honeycomb lattice. 



FIG. 1. Structure of a layer in (NBu4)2Mn 2 [Cu(opba)]3-6DMSO-H 2 



The partition function of Eq.(g) will now be treated by the standard techniques of high temperature expansion and 
Monte Carlo simulation to extract quantitative information on the system. 
In zero magnetic field, the partition function becomes 
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Z(T,Q)= I Vs \ Y[ 2cosh(i/3JS'||s l 

J \<ij> 

We get an effective ferromagnetic interaction between the classical spins : 

-fiH<&= ln[2cosh(i/3JS'||s i + s i ||)] 

<i,3> 

We can perform a high temperature expansion of the partition function of Eq. (|^) . By using the star graph technique 
P|, we have derived the expansion of laZ(T, 0) up to the 30th order in the variable K — ^/3JS. For the specific heat 
we obtain the following result : 
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According to Eq.(j^), the zero-field susceptibility, defined by \ — y qjj 2 
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can be expressed as : 
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where Wy = \t3JS\\ Sl + Sj \\, Wg = -±3JS(sf + s|), = tanh(W ij ) Wrf/Wij and 
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The high temperature expansion of x can be obtained in two i: 
function of Eq.(p|) both in powers of K and H and retaining the 



be obtained in two independent ways : first by expanding the partition 
w a^A ™ti™™n- coefficient of H 2 ; second, through an expansion of 
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the correlation functions which occur in Eq.(g|). Up to the 7th order, we apply both methods in order to validate our 
results. By using the second approach, more tractable at higher orders, we have obtained the following series, up to 
the order 11 : 
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In order to improve the range of validity of the expansions (Eqs.^jTj)) we performed a Pade approximant extrapo- 
lation toward low temperatures (large K values). For both series we get good stability of the Pade table. The results 
will be presented below. 

We performed a self-consistent check of our results by means of a Monte Carlo simulation of the effective classical 
model (Eq.(il)). The various observables can be expressed as ensemble averages with respect to the Boltzmann weight 
l/Z(T,0)e-P H ° t{ . For instance, if we define E = -k B TJ2 <hj> tanh(VF„) which is the energy of the quantum 
spins for a fixed configuration of the classical ones, we find that the internal energy is simply given by (E) H ff and 
the specific heat by 



C v = k B P 2 
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Similarly, from Eq.(g) the susceptibility can be expressed as a sample average. We have used the Wolf algorithm || 
adapted to our effective Boltzmann weight, on lattices of size increasing with K, up to 2 16 hexagons for K = 5. 

In Fig. 2 we have plotted the specific heat as a function of K. The data points correspond to the Monte Carlo 
simulation and the continuous line to the highest order diagonal Pade approximant of the high temperature expansion. 
The first remarkable feature is the well marked knee-hump variation of Cy as the temperature decreases. One can 
understand this effect in the following way. At very low temperature, the system is dominated by the effective 
ferromagnetic interaction between the classical spins. For the purely classical Heisenberg model, one expects a peak 
at low temperature in the specific heat corresponding to the crossover between the low temperature critical 
regime and the high temperature uncorrelated one. 
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FIG. 2. Specific heat Cv/NLks versus K — 1/2/3 JS : The Monte Carlo results are represented by squares, the solid line 
shows the Diagonal (7,7) Pade approximant corresponding to the high temperature series of Eq.(5) and the contribution of the 
quantum spins obtained from Eq.(9) is shown as a dashed line. The zero-temperature (K — ► 00) limit, Cv /Nhkn = 2/3, is 
indicated. 
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The hump observed in Fig. 2 near K — 4 corresponds to this effect. As T increases, the classical spin system goes 
rapidly to a disordered state, whereas the quantum spins remain locally coupled to their classical neighbours. The 
knee at K = 1 corresponds to this local antiferromagnetic order. This effect can be quantitatively confirmed by the 
following calculation. Assuming that the classical spins are completely random leads to (-E 2 ) = {E) 2 in Eq.(||) and 

C v = k B (<S>) disoidered = 8N L k B K 2 ^ —^—— dy (9) 

This function, displayed in Fig. 2 (dashed line), exhibits a peak exactely under the bump observed in the full C v . 

The other feature which emerges from Fig. 2 is the spectacular agreement of the series expansion with the Monte 
Carlo results up to K ~ 3 and the ability of this series to reproduce the double-bump structure. 

We obtain the same perfect agreement of both our methods for the magnetic susceptibility up to K = 3. In 
Fig. 3, we have plotted (solid line) the (6,6) Pade approximant of the series Eq.(^) as a function of temperature, in 
a range which corresponds to < K < 1.2. The data points correspond to the experimental results which will be 
discussed below. The shape of this curve can be explained as follows. Here again there are two physically different and 
competing correlation effects. The first one, due to the antiferromagnetic spin compensation, leads to a decrease of the 
susceptibility, more pronounced at lower temperature. The other contribution, due to the ferromagnetic correlation 
of the classical spins, induces a divergence at T — 0. 
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FIG. 3. Magnetic susceptibility times temperature yT in units of cm 3 : K : mol -1 versus T in Kelvin: The experimental data 
(squares) are taken from [9]; the solid line corresponds to the diagonal (6,6) Pade approximant based on the high temperature 
series Eq.(7) with J = 47.6 K, g 1 =2.0, g 2 =2.2. 



According to universality, we expect the critical behaviour at zero temperature to be described by the non-linear 
a- model. At low temperature the effective interaction Eq.(Q) reduces to 



■\JS ^2 II s * + s ill ~ const + g JS ^2 ® 
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which corresponds to the low temperature limit of the ordinary classical Heisenberg model on the honeycomb lattice 
with the ferromagnetic exchange constant J* = jJS. The long wave-length expansion leads us further to the total 
effective energy in the form of the non-linear cr-model Hamiltonian 
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(10) 



where n(r) is the three dimensional unit vector field on the two dimensional plane r. As an immediate consequence 
from the quadratic expansion we have the magnon contribution to the total energy per classical spin which is linear 
in T: 

E re -^JS + k B T 
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leading to the constant specific heat at low temperature Cy = ^N^ks- One can see in Fig. 2 that our Monte Carlo 
data are compatible with this value. 

As another consequence, we can borrow the low temperature behaviour of the magnetic susceptibility from the 
renormalisation group results [^|Jl2| for the non-linear cr-model of Eq. ( [Toj ) 

X = const • T 3 exp 



4tt JS 
k B T Ay/3 



const • T exp 



U/3 . 



(11) 
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FIG. 4. In (x# 3 ) versus K: The squares correspond to Monte Carlo results while the solid line shows the prediction of 
Eq.(ll) of the non-linear cr-model. 

In Fig. 4 we have plotted the low temperature Monte Carlo data for \n(K 3 x) as a function of K. The straight line 
corresponds to the behaviour of Eq.(|ll]). We obtain very good agreement, thereby confirming that at our lowest 
temperatures, we are well inside the expected universality region and the size of our simulated lattices is large enough. 

These results show that our series expansions constitute a reliable parametrisation of the specific heat and the 
magnetic susceptibility of the model up to K w 3. We can now apply this parametrisation to the experimental data 
H We have fitted the product T\ as a function of the temperature in the range 60/-T < T < 30QK with J, <?i and 
g2 as free parameters. The result is shown in Fig. 3, with the best-fit parameters given by J = 47. 6K, g\ = 2.0 and 
r/2 = 2.2. We observe excellent agreement with the experimental data, confirming that the high temperature magnetic 
properties of this compound are well described by our model. Furthermore, these values of the parameters are very 
close to those obtained for both Cu-Mn pairs |]l3| and chains [p|Jl4| with the same bridging network. The weakly 
pronounced minimum observed in Fig. 3 for T i=s 120K has its origin in the local antiferromagnetic ordering discussed 
above. 

The material exhibits a ferromagnetic transition at T c = 15K that the isotropic model cannot describe. To give an 
account of this nonzero-temperature critical behaviour the Hamiltonian must be generalized to include spin anisotropy 
and three-dimensional effects |l5| ], a modification which we are now investigating. However, it is amazing that the 
agreement between the experimental data and the model persists down to rather low temperatures (T w 20K), close 
to the measured T c . 
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